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Abstract 



We have performed a measurement of the i^ + ^7r + 7r°7 decay and have 
observed 2 x 10 4 events. The best fit to the decay spectrum gives a branching 
ratio for direct photon emission of (4.7 ± 0.8 ± 0.3) x 10~ 6 in the vr+ kinetic 
energy region of 55 to 90 MeV and requires no component due to interference 
with inner bremsstrahlung. 
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The radiative decay K + ^7i + TT 0, y (K^-y) in which the 7 is directly emitted (DE) is sen- 
sitive to important aspects of the low energy hadronic interactions of mesons Q. Although 
the DE component is difficult to observe because it is very small compared to the dom- 
inant inner bremsstrahlung (IB) associated with the K + — > 7r + 7r° (K^) decay, it may be 
isolated kinematically. The DE component consists of magnetic and electric transitions. 
At 0(p 4 ) in Chiral Perturbation Theory (ChPT), the magnetic amplitude consists of a re- 
ducible chirally anomalous amplitude Pj which is unambiguously determined by the Wess- 
Zumino-Witten functional || and of a direct anomalous amplitude [[|, arising from chirally 
covariant odd-parity octet operators. The direct anomalous amplitude is not calculable in 
a model-independent way but the size is expected to be comparable to that of the reducible 
amplitude or smaller. At 0(p 6 ) in ChPT, vector meson exchange (VME) contributions are 
also possible but, in contrast to the Kl — > 7r + 7r _ 7 decay 0, the VME contributions are 
less important; a weak deformation model |J predicts that the VME and the direct 0(p 6 ) 
contributions cancel [[|. There is no definite prediction from ChPT on the electric tran- 
sition amplitude, which depends on undetermined constants. Since the electric amplitude 
interferes (INT) with the IB amplitude, it may be distinguished from the magnetic ampli- 
tude, which does not. The experimental determination of the electric amplitude is of great 
interest, not only for ChPT, but for searches for possible non-Standard-Model effects like a 
CP- violating asymmetry between K + -^7r + 7r°7 and K~ ^7r~7r°7 decay widths predicted in 
supersymmetry J7|. 

Previous experiments |8|- [10|1 used decay-in-flight techniques to measure the branch- 
ing ratio in the tt + kinetic energy region 55 MeV < T + < 90 MeV with Abrams et al. || 
confirming the theoretical prediction for the IB branching ratio of 2.61 x 10~ 4 M. The 
current Particle Data Group average for the DE component of the branching ratio is 
(1.8±0.4)xl(r 5 
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which is five times larger than calculated from the reducible anomalous 
amplitude alone and consistent with the coherent combination of that amplitude and the 
direct anomalous amplitude suggested by a factorization model M . A decay spectrum shape 
analysis of Ref. || indicates the absence of electric contributions to the DE component but 
the existing data are not precise enough to rule out constructive interference with IB. 

The differential decay width for the decay K + (p) — > 7r + (p+)^ \po)l{q) is conveniently 
expressed in terms of a Lorentz invariant variable W 2 = (p ■ q){p+ ■ q) / (m 2 r+ m 2 K ) , where p, 
p + and q are 4- momenta of the K + , ir + and 7, and m n + and itlk are masses of tt + and K + , 
respectively. The differential width can be written in terms of the IB differential width as 
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where T + is the kinetic energy of the charged pion in the kaon rest frame, A is the decay 
amplitude for K + — > 7r + 7r°, and E and M are electric and magnetic DE amplitudes, respec- 
tively. The direct emission amplitudes E and M are independent of W and their effects 
are detectable only in the large W region, while the IB decay amplitude is proportional to 
1/W 2 , and thus dominates in the small W region. Experimental determination of these 
amplitudes requires high sensitivity to deviations of the measured spectrum from that of IB 
over a wide W region. 

The E787 detector ]12] at the Alternating Gradient Synchrotron (AGS) of Brookhaven 
National Laboratory (BNL) was used to study K n 2-y- Kaons of 790 MeV/c |13| were incident 
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on a solenoidal spectrometer with a 1.0 Tesla field at a rate of 7 x 10 6 per 1.6-s spill of 
the AGS. The spectrometer contains beam detectors, a scintillating- fiber target where the 
kaons stop and decay, a central cylindrical drift chamber, and a range stack (RS) of plastic 
scintillators with embedded straw chambers. Measurements of momentum, kinetic energy 
and range of charged decay products were performed. The output pulse-shapes of the RS 
counters were recorded providing good timing accuracy. A hermetic calorimeter system 
consisting of barrel (BL) and endcap (EC) detectors measured the positions and energies of 
photons and other particles from K + decays. 

The BL calorimeter was used to detect the three 7's in the decay. It consisted of 48 
azimuthal sectors and 4 radial layers of sandwiches of lead (1 mm) and plastic scintillator (5 
mm) sheets, 14 radiation lengths in depth covering a solid angle of about 3ir sr. The visible 
fraction of the energy deposited in this system was 30%. The Z position (along the beam 
axis) of the BL hits was obtained using the charge and timing information from phototubes 
on both ends of the 2-m long modules. 

The trigger for required a K + decay at rest, a charged track in a kinetic energy 
region between the endpoint of the K + — > 7r + 7r°7r° (K^) decay and the peak of the 
decay, at least three BL clusters, no extra particles in the RS, and no energy in the EC. The 
total number of kaon decays at rest available for the K^-y study was 1.8 x 10 11 , and a total 
of 1.1 x 10 7 events survived the trigger. 

The energy and direction of the three 7's from the decay were determined from 

reconstruction of the BL hits and the K + decay vertex position in the target. The 7r + 
momentum vector and the kinetic energy were measured with the target, the drift chamber, 
and the RS. Thirteen observables were available for the kinematic fit in this analysis. The 
relation between range measured in the RS and momentum measured in the drift chamber 
was examined to reject fi + backgrounds. Events were also rejected based on the presence of 
any activity other than three 7's and a 7r + . 

A kinematic fit was applied to the observables with the following six constraints: total 
energy-momentum conservation, consistency of the charged particle's energy and momentum 
with a pion hypothesis, and invariant mass of two 7's being equal to m^o. Since there are 
three possible pairings of the 7's to form the ir°, incorrect pairing can move IB events into 
the region of large W although IB events tend to populate the region of small W. 

To minimize the contamination arising from a wrong combination, the square of the 
matrix element of the IB decay was taken as a weight in the selection of the two 7's to be 
assigned to the 7r°, and the combination which maximized the product of the x 2 probability 
of the kinematic fit, Prob(x 2 ), and the IB weight was chosen among the three possible 
event topologies. Figure |I] shows the x 2 probability distribution of the fitted events. The 
distribution is consistent with measurement resolution throughout all values of W . Most of 
events at Prob(x 2 ) = are due to pion nuclear interactions and muon background, which 
do not satisfy the kinematic constraints. 

In the Monte Carlo simulation of the IB, DE and INT spectra, K n2 -y events were generated 
with the corresponding matrix elements in PJ and analyzed with the code used for real 
data. Thus, the effect of possible incorrect pairing of the two 7's from tt° on the Monte 
Carlo spectra was taken into account. The probability of incorrect pairing with this analysis 
was estimated to be 2.6% for IB, 20% for DE and 9% for INT. For IB, 2.3% of the events 
populate the region of W > 0.5, and only 6% of them were due to incorrect pairing. The 
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estimated W resolution was about 0.02. 

For the signal selection we required the \ 2 probability to be more than 10% and the fitted 
7r + momentum to be between 140 and 180 MeV/c; the latter condition on ir + momentum 
was imposed to remove the K n3 and backgrounds as well as K + — > 7r + 7r°7 decays outside 
the region 55 MeV < T + < 90 MeV. A total of 19836 events survived all selection cuts, eight 
times more than the previous experiments. Figure shows the 7r + momentum distribution 
of the K^-y events after applying all selection cuts. The 7r + momentum distribution is 
insensitive to the DE contribution, and thus it is reproduced by a Monte Carlo simulation 
of the IB contribution alone. 

Figure [|a shows the W projection of the signal events in the T + signal region. The 
dashed curve is a best fit to IB alone with a x 2 °f 61 for 7 degrees of freedom. The solid 
curve is the best fit to the sum of IB and DE and gives a x 2 of 7.9 for 6 degrees of freedom. 
Normalized to the IB spectrum the deviation of the measured spectrum at large W is shown 
in Figure ^|b. The DE component was measured to be (1.8 ± 0.3)% of the IB component 
in the region 55 MeV < T + < 90 MeV. The branching ratio for DE is determined to be 
B T+ (DE) = (4.7 ± 0.8) x 10~ 6 in the T + signal region by normalizing to the theoretical 
prediction for the IB branching ratio. The best fit to the sum of IB, DE and INT gives the 
INT component to be (—0.4 ± 1.6)% of the IB component with a \ 2 of 7.8 for 5 degrees of 
freedom. We conclude that no interference term is required for the fit to this spectrum. 

Several checks to verify the result were made. Spectral analyses were performed on the 
cos6' 7r + 7 distribution, where 9 n + y is the opening angle between a 7r + and a radiated 7, and 
on the energy distribution of the radiated 7. The resulting branching ratios obtained for DE 
are consistent with the W fitting result to within 3%. This difference is used as an estimator 
of the systematic error due to the fitting method. Selection of 7 pairs without IB weighting 
reduced the branching ratio by 6% while doubling the statistical uncertainty. Systematic 
uncertainties due to calibration and resolution of the energy and position measurements 
were estimated from the observed DE branching ratio variation as each corresponding Monte 
Carlo parameter was varied over an expected range. The uncertainty due to limited Monte 
Carlo statistics was estimated by performing the analysis with small sub-sets of the Monte 
Carlo data. All systematic uncertainties in the DE branching ratio were added in quadrature 
to obtain an overall relative error of ±6% (±0.3 x 10 -6 ). Table | gives a summary of the 
systematic uncertainties. 

Two other possible background sources in the large W region are K + — > 7r°/x + i/ M decay 
with a radiated or accidental 7, and K nS decay with a charged pion decay- in-flight. Both 
processes contain a charged muon in the final state. The estimated number of the muon 
background events was 103 in the whole W region (7 out of a total of 597 with W > 
0.5), which was obtained by selecting events with a muon-like range-momentum relation. 
Imposing additional pion identification criteria based on the presence of a 7r + — ► decay 
signal in the stopping counter [jnj did not change the W spectrum. The insensitivity of the 
final result to the variation of photon veto cuts also confirmed that these backgrounds are 
negligible. 

Assuming a pure magnetic transition for the DE amplitude, the result for B T+ {DE) gives 
\M\ = (2.1 ±0.2) x 10 -7 for the dimensionless magnetic amplitude in Eq. (1). The reducible 
anomalous amplitude | — eGgm 3 K /27i 2 f\ in Ref. is \M\ = 1.8 x 10~ 7 with standard 0{p 2 ) 
ChPT coupling constants, / = 93 MeV and |G 8 | = 9 x 10~ 6 GeV -2 . This result supports 
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the hypothesis that the dominant contribution to DE is due to the reducible anomalous 
amplitude and the other magnetic contributions are small or cancel. Similarly, the null 
measurement of the INT component gives a 67 % confidence-level interval of —2.6 x 10~ 8 < 
Re(E) < 1.6 x 1(T 8 for the electric amplitude. 

In conclusion, we have performed a new measurement of the K + — > 7r + 7r°7 decay with sig- 
nificantly higher statistics and improved kinematic constraints, using the E787 detector and 
K + decays at rest. The branching ratio for direct emission is 4.7±0.8(stat.)±0.3(sys.) x 10~ 6 . 
The decay spectrum indicates that the direct emission amplitude is consistent with being 
due to purely magnetic contributions comparable to the reducible anomalous contribution. 
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TABLES 



Estimated systematic uncertainties 


(%) 


Fitting Method 


±3 


7r + Momentum Calibration 


+3 

— 1 


7 Position Calibration 


+3 
-0 


7 Position Resolution 


+3 

-0 


7 Energy Calibration 


±0.3 


Monte Carlo Statistics 


±2 


7 Interaction in Material 


+0 
-4 


Total systematic error 


±6 



TABLE I. 

A summary of estimated systematic uncertainties. All uncertainties are added in quadrature 
to get the total. 
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FIG. 1. Distribution of the x 2 probability for the fitted events. Events with Prob(x 2 ) > 0.1 
are accepted as due to K + ^7r + 7r°7 decay, (a) All events, (b) events in a region of W > 0.5. 
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FIG. 2. Momentum distribution of observed tt + events (dots) and Monte Carlo simulation 
with IB alone (solid histogram). A \ 2 of the fit is 18 for 15 degrees of freedom. 
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FIG. 3. a) W spectrum of the observed events and best fits to IB+DE (solid curve) and IB 
alone (dashed curve); b) W spectrum normalized to the IB spectrum. 
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